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BMP SPEUCOJIffiE CMCULATIOif OP EIEGPRON MOMMTUM 
DISTEIBUTIOIIS IN SOl'IB TEMSITION MEIiiS 

By 

DILIP GOVIND KMHBBl 
Itepartment of Physics 

Indian Institute of Technology, Eanpur, India 

May 1976 

In recent years many detailed studies of the electron 
momentum distribution (EMD) in solids have been ma^de by experit 
mental techniques employing gamma— ray Compton scattering and 
positron angular correlation measurements. A p 2 X>per comparison 
between theory and experiment is best made through single 
crystal profiles and their anisotropies. Simple theoretical 
models are not useful for calculating EMD in crystalline solids 
and, therefore, energy band calculations beccme necessaiy. Such 
studies are of particular interest for transition metals in 
view of their interesting electronic structure. The present 
thesis describes a band structure calculation of the EMD (for 
Compton scattering) in Pd and Ag employing Hubbard-Mijnarends 
method. The emphasis in the present work is to examine the 
role played by band structure in determining EMD, A systematic 
study of the EMD in the 3d- and 4d- transition metals is also 
made. In the end a computational procedure to calculate 
directional anisotropies tTcm theoretical is developed. 
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In Chapter I the techniques of Compton profile (CP) and 
angular correlation of positron aimihilation radiation (^PAR) 
are introduced eind am outline of relevant theory and experimental 
methodology is given. Yarious problems connected with the 
determination of IS-ID in metals are discussed. This is followed 
by the compairison of these two techniques. 

A renormalize<i free atom (RFA) model for the calculation 
of CP and ACPAR is discussed in Chapter II. A calculation of 
GP and ACPAR for Ni and Cu is presented and the results are 
discussed. It is found that the RFA model gives a good overall 
description of the GP data for Ni and Cu and for Ni it supports 
an electronic configuration closer to (3d)^(4s) . As far as the 
ACPAR data are concerned the RFA model does not give an equally 
satisfactory description. The need for a bemd structure calcu- 
lation of EMD is emphasized . 

Chapter III presents the formalism of Hubbard -Mi jnar ends 
band structure method of calculating IHD in transition metals. 
Necessary formulas and computational details are given and the 
limitations of this method are pointed out. 

A calculation of the band structure and IMD (for CoH^jton 
scattering) for pd and Ag obtained by en^loying above method 
is presented in Chapter IV. The results for the energy bands 
are sho^m to be in satisfactory agreement with other ri^orted 
calculations. Band-by-bamd contribution of the IMD is shown 
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along the three ^ and A) symmetry directions and discussed 
in the light of respective band structure. In order to gain 
more insight, constant EMD contours (isodenses) in the (lOO) 
and (110) planes for total H4D as well as for separate bands 
are plotted. These contoux’s ar^e analyzed and compared, bringii^ 
out the role of d bands, Fermi surface etc. in the EMD. Results 
for the iso elec tronic pairs Ni and Pd, and Cu and Ag are coiapared. 

In Chapter V, a systematic study of EMD as one goes along 
the 3<i- and 4d- transition series is made. Ihe effects on EMD 
of the different properties of the d-band, e.g, bandwidth, s_d 
hybridization etc, are discussed. A ccHnparative study of EMD 
in 3d- and 4d- transition series is also made. 

In Chapter VI, a method of calculating single csrystal CP 
(or ACPAR) profiles from theoretical EMD enploying an expansion 
in termis of Kubic harmonics is developed. For Cu, the method is 
shown to give qualitative agreement with CP and ACPAR results. 

The anisotropies of Pd and Ag are predicted and discussed. 

Various aspects of this method and its limitations etre discussed. 

These studies have enabled the author to gain valuable 
ij^ight into the band structure effects in the EMD of transition 
metals. It is hoped that these studies have contributed to our 
unders tanding of the inle played by band structure in determining 
EMD, particularly in the transition metals. 
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1.1 Int-roduc-fcion : 

During t±ie last few years the electron momentum distri- 
butions in atoms, molectiles and solids have been a subject cf 
Considerable theoretical and e 3 q>erimental interest, Ihe first 
direct experimental evidence for the dynamical nature of atxaic 
electrons came in 1928 from the observed broadening of the 
inel astically scattered X-rays from atoms, DtMond m pointed 
out the importance of sxich experiments for stxxlying the electron 
momentim distributions (EMD) and ascribed the broadening of the 
Compton-scattered radiation lineshape (the Doppler broadening) 
to the linear nv^nentm of the ^ectrons. 

The next three decades saw no significant development 
of this experimental technique, the main reason being the ron- 
availability of highly intense X— ray sources. It was in 1965 
that Cooper, Leake and Weiss ( 2} demonstrated the feasibility 
of such experiments for studying EMD of metals like Li, Studies 
of other systems followed in the next few years. However, the 
most significant development occurred in 1971 when gaamna-ray , 
Ccxnp ton scattering "^s introduced f 3 |, During the last few 
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years the growth of this technique has been rapid^ and today 
it is anerging as a powerful tool for investigating IKD in a 
wide variety of materials [4-7]. The interesting history of 
these developments is outlined in a comprehensive review by 
Cooper [ 5], # 

During the period the technique of Compton profile lay 
dormant, another technique to stxady EMD was established and 
widely applied using a different microprobe, viz. positron [ 8-10] . 
The feasibility ©f measuring angular correlation of positron 
annihilation radiation (ACPAR) was demonstrated [ill in 1950* s. 

In the last decade many acciarate measurements on several 
substances employing ACPAR technique have been reported 
indicating the potential this technique for investigating 
EMD in solids, in particular for studying the Fermi surface (FS) 
in metals [ 10 L 

In both these experiments the momentxm of the electron, 
rather than its energy, plays a more fundamental role, Ihe 
momentum of the electron is directly related to the ground 
state wavefunction of the system and therefore the knowledge 
and analysis of the EMD should give us interesting information 
about the ground state properties of the electrons, thus providing 
a test of the accuracy of the electronic wavef unctions. Parti- 
cularly for metals one expects to obtain valuable information 
about the Fermi surface* One must, however, iK>te that it is 
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not the EMD itself ■which is measiired in -the alDove experiments 
b'ut a line or surface integral of the SMD in the momentum space. 
The positron and -the plxs-bon interact wi-th the electrons in 
different ways and therefore the Comptaon scattering and the 
positron annihilation techniques differ in -Cerms of "the experi- 
mental acc’uracy offered^ -theoretical -understanding gained and 
to some extent the type of information ob-tained [9], However, we 
believe -that -these techniques are complenentary to each o-ther 
and a greater xjmder standing about the electronic states can be 
gained by studying a particular system anploying both these 
techniques. We shall now discuss both -these techniques in some 
detail, restricting o-ur discussion mainly "to those aspects which 
are directly related to EMD, 


1*2 The Compton Profile ; 

1 * 2. 1 The Basic Process : 

The inelastic scattering of axphoton by elec-trons is 
known as Compton scattering. It is well-kno-wn that for a 
stationary elec-tron the shift (AX) in the wave length of a 
photon, Comp-bon-scattered thro-ugh an angle ♦ , (Fig. 1,1 ) is 
given by. 



^ sin^ ^ 
me 2 


( 1 . 1 ) 


However, -the above rela-fcion is modified in the case of a free 
el ee-tron having a finite linear momentvHn, So -that -the 





1.2: Schematic diagram of the Comptnn profile. 
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Compton-scattered line is Doppler— broadened. The energy momentim 
conservation laws [5 } give a modified escpression. 


AX = 


■ 2 f 

me 2 ” me 


2h 


Sin 


( 1 . 2 ) 


where ^ is the wavelength of the incident photons, p^ is the 
z— component of the electiron momentum p along the scattering 
vector (which is taken to be z-axisO, It is seen from (1,2) 
that the intensity distribution, I ( x) , in the scattered radia- 
tion lineshape will be determined by the z component of the 
electron momentum. Such a broadened lineshape or the plot of 
l(x) versus X (Fig. 1.2), is generally krtown as the Compton 
profile (CP), We shall now briefly discuss the theory of CP* 
in the impulse approximation (lA), 


1.2,2 Impulse Approximation ; 

Consider the inelastic scattering of a photon from a 
system containing electrons. Let k^ and be the wave vector 
and frequency of the incoming photon while k 2 and represent 
Corresponding quantities for outgoing photons, A schematic 
diagram of the prodess is already shown in Fig, 1.1. Ihe 
Haniltonian for an electron in an electromagnetic field is giv^ 


H 


•*■2 2 
+ V{r) + 

2"^ 2mc^ 


e T 

me 


(1.3) 


where for simplicity we have considered only a one— ^ectron 



6 


system# In (I^B) tlie terms in A and p#A are to be treated 

as perturbations. Since the scattering is weak the cross 

section for th.e process can be calculated in the lowest order 

Bor^i approx in ation, I^e terms responsible for the scattering 
, ~*2 

are the A term in the first order and p,A term in the second 
order. It can be slno-wn [12] that for energy transfers large 
compared to the bincdirg energy of the electrons the contribu- 
tion from the p,A tern is small and can be neglected, Ihe 
crDss~sect±oTi for tbe Compton scattering (ti = c = l) is then 
given by. 


dfld 

where (da/dSil 


_ » Thojisom cress-section 
Th - , 

r e ■ \2 V 


ilc. 


{ i>| S (E^ - 


— ti) ) 


( 1 . 4 ) 


|l> and |f > are respectively the initial and final states of 
the electron^ •wbil.e Ic = — k 2 snd cu = are respectively 

the momentm and energy txansferred to the electron. The ® 
-function represents energy conservation and and ^ denote 
the polarization" ■vectors for the initial and final states of 
the photon. 

We can rewrrite (1,4) as 


d g 
d£ld® 




( 1 . 5 ) 



/ 


where, 

sat, to) =l {<f { ji>j^ 5(S£ - 

f 

- 4 * * 

and S(]c, ) con-tains all the information about the state of the 

electron and is Icnown as the dynsmical scattering fac-tor. It is 
related to the EMD and is thus of interest -to us. 


Now we assime that -the energy transferred to the electron 
during the collision is large coirpared to the binding energy of 
the electron. In other words, we assume that the time of the 
interaction is short enough to consider the particle unbound 
during the scattering (i.e. the particle sees a constant poten- 
tial during the interaction). This is the essence of inpulse 
approximation. Under this approximation (IA), S(k, “) can be 
written as. 


where 

that 


S^Clc,«i)) = J dp*" p(p)- 5(w- ^ - ~^) {l*7a) 

= J- (q) (1.7b) 

}c k 

p(p) is -the EMD and J*(q) is the Cotapton profile (CP) such 

k . 


P(5) = 1 X (p) P = -^j drf' Xl.S) 

( 2 ^) 

and 

J.Cq) = f dp p(p) 5(q -p.k) (1.9) 

k 
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where q = 7 ^ 7 - — ^ lk| , and the quantity q is thus the measure 

Ikl 

of the electron moment'cm in the direction of momentun transfer. 
If z—axis is chosen along k, the scattering vector, the e>qDres- 
sion for CP becomes, 

- 1-00 

J*(q) = IJ dp^ dpy [ P(^]p (1.10> 

z 

— CO 

Equation (l,7a) clearly shows that the scattered photon 

,2 

is shifted in freauency both by a momentum transfer term 

i* 2 m 

and a Doppler shift term (-^^) , It is seen from (1.9) that the 
CP is simply a projection of the EMD, p on the scattering 

vector. The above form (eq, 1,8) of P(p) can be extended to a- 
many-electron system, by writing, 

. occ ^ 

P(p) = I % ( 1 . 11 ) 

where f^ is the occupancy of the state i. 

The simple and straightforT.'rard relation (l.lO) between 
Etm and CP rests on the validity of the inpulse approximation 
which is, therefore, very important for the determination of 

1,3 Validity of the Impulse Approximation : 

The validity of the lA is studied in scane detail by 
Eisenberger and Platzman f 12 3 , and by Currat et al, El3 ] , 

By considering Is hydrog^ic systems it is sh>wn [l23 that 
the correction to lA is of the order of (Eg/Sj^) where Eg 
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is the binding energy and is the recoil energy of the 
electron. For a typical X— ray energy { “20 keV) the correction 
to the lA is quite small for the outer electrons because 
Eg “ 10-100 eV. A calciolation for Li (is), C(ls) and Al(2s/-2p) 
employing CKCited state one— electron continuum "wavefunctions 
showed [13 ] a good agreoment with the results onploying lA. 
Exanination of the validity of lA and comparison with experi- 
mental resxiLts indicate that lA adequately describes CP through 
(l,lo). Our calciolations described in Chapter 6 are made within 
the framework of I A. 

1,4 Experimental Determination of CP t 

An ecxperiment to determine the CP consists of a measure- 
ment of the energy spectrum of photons Compton— scattered at a 
fixed angle 0 (Fig. 1,1 and 1,2). Earlier measurements of CP 
were carried out employing X-rays, from either Ag K (22,2 keV) 

ct 

or Mo (17,4 ke^) radiation. Experimental arranganent for 
X-ray Compton scattering is described by Cooper [5], In 1911 , 
Felsteiner et al. [3] introduced 59,54 — keV gamma rays from 
together with a Ge'(Li) spectrometer to study CP, Since 
then the 159.0 - keV garrma rays from Te [ 4] as well as the 
412 — keV gamma rays from Au [l< } have also been employed. 

The details of the gamma ray Con^ton scattering experiment and its 
methodology are described in the literatrsre l 4, 7,9] and today 



almost all the CP measurements are carried out with ganana rays. 
There are some definite advantages offered by the gamma ray mea- 
surenents as compared to the earlier X-ray technique [4] , 

These are: 

i) The gamma ray sources give monochromatic radiation as I 

compared to the X—ray sources which usually emit a do^iblet 

Ce.g, K ) riding on a continuous backgroxond of Broum'oirrahiamg ^ 

H “11 : 

ii) .The entire energy spectrxm of the scattered gamma rays io 
accumulated almost simultaneously, unlike the point— by— point 
wavelength scan in the case of X— rays. This resiilts in the 

i 

reduction of coxjnting periods. I 

iii) The higher energy of gairana rays improve the validity of | 

impulse approximation* allowing the interpretation throxigh (lilO), { 

Secondly the ratio of Compton to plxjtoelectric .cross— section goes s 

3 

as (Ey/Z^ , The use of ganma rays give a better Compton to 1 

photoelectric ratio even for high Z and make it possible to [ 

study elements with high Z, p 

Data from escperimentally measured CP have to be | 

t 

processed for effects like detector resolution, multiple scatter- ! 
ing, background correction etc. Systematic studies of these 
effects during the last few years [7] have enabled us to carjy 
out such a data analysis with proper understanding* | 

■ I: 

. ■ i 

We shall now briefly discuss another in^artant technique j 

' !■ 

Cviz, ACPAR) used for studying iKD, | 



1.5 Positron Anqtalar Correlation Stxidjes: 


Tlie use of positron annihilation process for investiga- 
ting solid state properties is novj well established [S], The 

low energy ( = 200 keV) positrons emitted from a suitable 

22 

radioactive source (e.g, Na) penetrate into the matter/ 
thermalize in short time ( “ lO”'^'^ sec) and then einnihilate 
with electrons. The annihilation characteristics depend 
almost entirely T: 5 )on the state of the positron and electron 
just prior to the annihilation. The information about the 
annihilation process is transmitted to the observer via the 
annihilation ganroa rays (through parameters like energy/ 
momentim and lifetime) without appreciable attenuation or 
scattering. These feat^ares make the positron annihilation 
technique attractive for stvidying the electronic properties 
of matter. 

The experimental techniques involving positron annihi- 
lation can be divided into following methods, 

i) Study of two-photon angular correlation, 

ii) Study of Doppler-broadening of annihilation 

radiation lineshapes, 

iii) Study of positron lifetimes/ 

iv) Study of 2Y/3 Y yield. 

The first two methods are directly relatai to the 
study of IMD, In the following we outline the first method. 
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Consider the process in which the positron annihilates 
with an electron having antiparallel spin giving rise to two 
gasmia rays. The law of conservation of momentm demands that 
the emitted photons be collinear in the centre-of-^nass systen. 

In the laboratory frame of reference, the annihilating 
positron- electron pair possesses finite momentim and the 
annihilating photons, therefore, deviates from collinearity. 

This deviation can be observed as an angle e^, = p^/mc which 
is usually of the order of a fev? milliradians (l a.u. =7.3 
milliradians). Since the positron is thermalised at the time 
of annihilation its momentim can be neglected. It is then 
mainly the momentun of the electron which gives rise to a 
deviation from collinearity. 

Considering the initial state of the systaii as a product 
of the ground state electron and positron wavef unctions, the I 

two-photon momenta density, P(p‘)/ is given [ s] by i 

2 ■ 

p(p*) = — f (?) (?) dr j ^ (1,12) 

where P(p*) dp' is the probability of annihilation with the [ 

3nission of a photon pair having a momenttm between p* and | 

p* + dp*, ^(r) and f_(r) are the ground state wavef unctions I 
for the positron and electron respectively, and = e /fee i 

is the classical electron radius. For the case of 'tiaermalised | 
positron, PCp* ) gives a good approximate description of the I 
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HMD itself. However the positron being positively charged 
tends to keep away from the nucleus and hence annihilates 
preferentially with the outer (or conduction) electanons. The 
expression (1,12) is identical (barring some constants) with 
(1,8) for the case of constant positron waveftmction. 

Unlike CP studies, experimental determination of 
angular correlation of positron annihilation radiation 
(ACPAR) employs a variety of geometries fS/lo] , In the 
so-called loQg— slit geometry one measures the coincidence 
counting rate N(p'), which is related to P(p*') through the 
following relation, | 

= // P^P'^ <SPx ^Py (1.13) 

Other geometries (e.g, point slit or mifLti -counter geometry) j 

s 

on the other hand measure a coincidence counting rate, i 

N (p^, = /^P(p’) dpx (1.14) I 

i 

An account of how different gecxnetries have been used I 

for the study of Fermi surface (FS) or EMD is given in 

literature [3,10]. | 

< I, 

From the above discussion of the two methods it is i 

' ■ ■ , ' . t' 

clear that the mcxHentrsn density, P (p) or p(p'), plays a central 

■ . ■ ■ , . 

role in the interpretation of the observed profiles (CP or ACPAR) ^ 
In the next section we shall briefly examine the problem of 
calculating IHD, 



1.6. GalciJlafclon of Momentun Density ; 


1.6.1 The Electxon Wavef\inct.ion i 

Computation of realistic electron wavef unctions in a 
metal is a major important problem in physics. Earlier 
theoretical calculations of EKD employed either a free electron 
or free-atom wavef unctions. Although a nxnber of band struc- 
ture methods have been available all these years, it is only 
recently that these methods have been extended to derive 
momentum wavef \anctions and momentum densities. A nuriber of 
calculations using APW, KKR, OP'-'J etc, are now available, 

1.6.2 The Positron Wavef ixnction i 

A knovrledge of positron wavef unction, ^(r), is indis- 
pensable for any meaningful interpretation of ACPAR data. 

The methods [8] commonly used for calcTolating ^^(r) are 
a) Wigner-Seitz method or b) plane-wave expansion method. 
Alternatively one can employ any of the standard band structure 
methods and calculate (r) for the lowest (ic = O) state. 

-r "T 

An interesting method developed by Kubica and Stott [ 15 ] employs 
the concept of pseudopotential for the positron and combines 
the merits of Wigner-Ceitz and plane wave expansion method. 

1.6.3 The PositTOn-Electron Correlations : 

Many authors C 16, 17 I have studied the effect of positron- 
electron correlations on annihilation rates and AGPiR curves 
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OTploying many-body theories. Such studies have indicated 
that the e"*" - e” correlations have a pronounced effect on the 
annihilation rates which depend on the effective electron 
density at the positron site. However, their effect on the 
momentim distribution c\irvo is less marked and, in particular, 
the position of the discontinuity in the moraentisn distribution 
(at p = pp) shows no shift due to these correlations [16]. 

This conclusion is especially inportant for the use of ACPAR 
curves for PS studies. An analogous result is obtained by 
Luttinger [igTin the case of interacting electron gas, 

1,6.4 The Electron-Electron Correlations : 

The electron-electron interaction can play a signifi- 
cant role in the interpretation of CP data from simple (e.g. 
Li, Na, Al) metals. The properties of the interacting 
electron gas are studied in some detail [19] and their effect 
on EMD has also been examined [20] . In our calculation of 
in transition metals, v;e have neglected the electron- 
electron correlation effects, apart from those entering via 
the one-electron potential.- 

In the following we shall briefly outline the form of 
EMD obtained in different cases. 
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1.6.5 The Free Electron Gas t • 

In the case of non-interacting electron gas the 
momentim density has an isotropic and constant value i 5 >to 
Fermi momentim Pp vdiere it suddenly drops to zero (Fig. 1.3a). 
In this case the CP or ACPAR is given by the relation. 


IT(p^) « 2ir (p^ _ p^) 


if < Pp 


= O 


if P 2 > Pp 


(1.15) 


v/hich has parabolic behaviour-* However, in the case of 
interacting electron gas electron— electron correlations 
push a part of the occijpied states below Pp -to the higher 
momenta values. As a consequence the sharp disco ntintiity 
smoothens slightly and a tail is developed in the momentm 
density (Fig. 1,3b), 

1,6,6 The Free Atom; 


The mom en tun distribution due to atomic s, p and d 
electrons (for the same shell) is shovm in Pig. 1,4. It 
is seen that only the s electrons can contribute to iMD at 
p = O. Hie reciprocality of the r - and p - space is also 
quite clear from Fig, 1,4. , The wavef unctions which are more 
spread in r - space are rrvore peaked towards the low-momentiro 
region in the p— space. The momentum distribution in this 
case, is isotropic and is given by. 
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Momentum >(a.u.) 
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where sore atomic wavef unctions. Usually the ’core* 

contribution to the El-iD in metals is calculated by using (1.16), 

1.6,7 The Conduction Electrons: 


In the case of metals, the conduction electrons, as 
v/ell as the positron states are Bloch states, Ihe wavef unction 
of an electron in a state 1c and band j can be written as. 





(r) 



(1.17) 


with u^ (r) = i (k) e 

k, j G G 


iG. r 


(1.18) 


vrhere G-^ are the momentm wavef unctions and G’s are the 
G 

reciprocal lattice vectors. In a similar fashion the plane 

wave expansion of positron wavef unction (r) would give 

k=o 




(r) 


k=o 


^_^(r) = 


I 


p 


iG’.r 


(1.19) 


G' G' 


The momentum densities for the two cases are then 
given by, 

p(p)« S p (^ 

KJ 

« ^ I I f - (f -t) 

Cj g ^ 


(l.,20a) 
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and 


p(p') , J- l\ cl (k) P 

k, j G- G* 


where 


(1?) 


I (k) 

G G G-G' 


(l.20b) 


(1.20c) 


In(l,20) the first sxjnmation indicates that all the occupied 
states below the Fermi level should be considered. The second 
summation^ together with the ® -function shows that electrons 
in state k will contribute not only at momenta p = k (normal 
or N— process) but also at momenta ^ = i? + (? (UmlcLapp or U— process) 
The contributions to P (^) from the states in completely filled 
zone will be continuous but that from partially filled zones 
will be non-zero only at p such that lies inside the 

FS. The momentum density then exhibits discontinuities not 
only at p = p^, but at p = p^, + G also/ its magnitude being 
weighted by a factor G^(1c) , Fig. l.Sc^ 1.3d, l,3e, and l,3f 
shows the behaviour of P(p) for some interesting cases. In 
this manner the effects due to i) Fermi surface ii) momentum 
wavef unctions and iii) U-processes combine to build the FMD, 

P (jx) and consequently the CP, J(q) or ACPAR, K(p ), curves. 

The inner ‘core' states being nearer nucleus give an 
appreciable contribution in the high momentum region (so-called 
high rocffnentum components, HMC), In the case of positron 
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annihilation, the overlap of positro -electron wavef iinction is 
appreciable at large r-values and the tMC arise mostly out of 
the U-processes. For Compton scattering, however, the I#IC 
derive their contributions from the core states as w^l as 
from the U-processes, 

Unless a pseudopotential method is used, one does not 
calculate the actual electron states in solids using plane 
wave expansion Cl. 17). Although the details of p(p) may 
change depending on the method used, essential foimi of p(p) 
remains same. 

Details of a band structure method for the calculation 
of fJ'lD will be discussed in Chapter 3. 

1.7 Comparison of CP and ACPAR Techniques ; 

A. comparison of the CP and ACPAR methods [ 9] is given in 
Table 1,1. From experimental point of view, the choice of 

various geometries and superior momcntu.i resolution are the 

i' 

main advantages of ACIAR method. liowcver, for theoretical 
interpretation, it is more straightforward to calculate CP 
which do not require a knowledge of positron wavef unctions. 

The CP and ACPAR methods can be complimentary?- -to each o-ther 
and it is desirable to experimentally stiady the same systans 
by both the techniques as well as to -theoretically calculate 
the two momen-tum densities employing s^e method and same 
electron wavef -unctions. 



22 


Table 1.1 

Comparison between positron annihilation and Compton 

profile methods 


Featxire 


Angular correlation of 
positron annihilation 
radiation with LS 
geometry 


Compton profile 


Probe 


Positron 


Photon 


Quantity 

measured 


Angular correlation of 
the annihilation photon 
pair lUp^) 


Information 
provided and 
assumptions 
of validity 


^(p^) = ^Px ^Py 

(i) Positron is therma- 
lized before 
annihilation; 

( i i ) po s i tro n-el ec tro n 
correlation is 
negligible 


Typical momentum O, C7 at. 
resolution # , 


units 


Counting 

efficiency 


Relatively lo^v'/ typical 
counting periods of 
1 week per (hkl) 
direction 


Quality of the Must be free from 
samples (single defects; surface purity 
crystals) very essential 


Data corrections Straightforv/ard 


Intensity distribution 
of the energy spectrum 
of the scattered photon 

JCp^) 

dPy 

impulse approximation 


0,4 at, xinits Pf 


Higher if -rays and 
Ge(Li) detectors are 
used; ti’^pical counting 
periods of 2 to 3 days 
per (hkl) direction 

Optimism thickness 
necessajry; too thin 
samples vsould give 
poor y i el d Sv whil e 
too thick sanples 
will cause multiple 
scattering effects 

Complicated. 
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1.7,1 Comparison Between Hieory and Experiment; 

As far as EMD is concerned two approaches are available 
for a comparison betv/een theory and experiment. First method 
consists of reconstructing the EMD, P(^)/ from CP (or ACPAR) 
curves measiired for a nunber of different orientations of a 
single crystal specimen, and comparing this reconstructed ES-tD 
with theoretical EMD. Alternatively, one can theoretically 
calculate a ‘model ' EMD- and obtain single crystal CP (or ACPAR) 

“ '^h*lc*l* (1.10). 

Although the former ‘unfolding* method, employing an expansion 

of in terms of lattice harmonies, as suggested by 

Mijnarends [ 2j Tias been tried on several metals, it has ix>t 
yet proved to be a unique metlxsd owing to the problems inherent 
in any ‘unfolding* or ‘deconvolution’. Perhaps a better 
strategy is to combine the above two methods. 

In view of this v;s have extended (Chapter 6) our 
calculations of E21D to predict directional CP curves and 
their anisotropies, 

1,8 Motivation Behind Our Present Work ; 

In this section we wish to lay down the motivation 
behind our present work. 

It is seen from above that for any meaningful comparison 
between theoretical and experimental CP or ACPAR curves, a 


curves J (q) or the anisotropies 
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calculation of the EI-ID is essential, Arong metals/ the EJID 

of transition metals is one of the most interesting because of 

their interesting electronic structure [22,23] . The best 
» 

approach for calcvilating the EJiD of transition metals is the 
band structxore method so that the effects on EMD due bo Fermi 
surface, the d bands, s— d hybridization etc. can be properly 
examined. The B4D in 4d transition metals hitherto appears to 
be an unexplored area except for the -worlc on Kb [24] and 
Mo [25. X, 

We foTJnd it interesting to calculate the EMD of Pd and 
Ag for the following reasons. In the first place the EMD of 
the 3d metals Ni and Cu, iso electronic with the 4d metals Pd 
and Ag have been well stiKlied theoretically [26] and experi- 
mentally [ 27 X , It should, therefore, be interesting to study 
the EEiD in Pd and Ag and compare it with Ni and Cu (all fee 
metals) to examine the effects of band structure, Fermi surface 
etc. Again, the band struct'ure and electronic properties of 
Pd are among the most fascinating . among the 4d metals [28}, 

In view of the fact that there are not many band structure 
calculations for the 4d metals, we thought it better to employ 
the fast but approximate method of Hubbard and Mijnarends [2^ 
rather than use a more sophisticated band structirre method, 
Calc\jlations of EMD in Ni and Cu employing Hubbard-^'Ii j narends 
method are already available l26] so a direct oompeirison of 
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these results is also possible. It is true that noK experi- 
mental data on Compton profiles of Pd and Ag is yet available 
for a comparison between theory and experiment but it is 
expected that such stxadies on heav\^ metals anploying 412— IceV 
gamma rays v/ould soon be made [14 ]. Until then our studies 
provide computer- simulated Compton scattering studies of Pd 
and Ag and can perhaps provide suggestions for designing and 
analyzing an actual experiment. In the absence of any experi- 
mental data^ the emphasis in our approach in the present work 
has been to study the role of band structure and Fermi sixrface 
on EMD, In other wordS/ we feel that EMD is yet another ground 
state property (e.g. electronic specific heat, susceptibility ^ 

etc, ) which can be predicted by a band structure calculation and | 

compared with experiment. In the past such calculations in ' 

I 

momentum space have been ignored and we feel that our calculations 
provide a systematic review of the results of HS4D in transition 
metals. To this end, we have discussed in Chapter 5, systeraatics ; 

of the SMD in 3d metals. Again a comparison between the ESdD | 

for Compton scattering, p(p), and the two— photon momentum dis- I 

, j 

tributions, p (p*), is very instructive and interesting. The I 

l: 

angular correlation of positron annihilation in Pd, Ag and 

Pd-Ag alloys are being studied theoretically and ^cperimen tally [29] 

' 1 

and v;e feel that our EMD calculations can prove interesting and I 
complementary data for such a ccanparison. i 
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Lastly we felt that a comparison of theory and 
experiment in this area can be made rwDre sensitively through 
the anisotropies “ '^100 than through single 

crystal profiles, ■^hkl- We have therefore extended our 
EKD calc^fLations in Chapter 6 to predict the anisotropies 
in Pd and Ag, 
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CkapteJt 2 


ELHmm lAOMEhTTW PJSmiBUTJOU W HTCKEL AW 
COPPER EhiPWyJUG REMORMALUEP-EREE^ATOU mOEL* 


2,1 Intxoductlon : 

In this Chapter we present an application of the 
renormalised— free-atom (R.FA) model to the calculation of 
HMD, The RFA model was developed by Hodges et al. [ij who 
applied it to calciilate some gmund state properties of flae 
transition metals, Berggren [2] employed RFA wavef unctions 
to calcxfLate the SMD for Compton scattering in Ti [ 3 ] , 

V [2 3 / a-nd Pe [ 4 ]and showed that the gross features of 
the CP data could be described by this computationally simpl 
model. We have extended this method to the calculation of 
EMD for positron annihilation and have calculated the CP 
and ACEAR curves in Ni and Cu, and compared them with the 
experiment. 

The motivation for the present worlc was as follows: 

1 ) To test whether the RFA rtodel is valid equally for 

describing positron annihilation and Compton scattering 
processes, • • 

* WDrk described in this Chapter is already px±ilished as 
a paper: D,G, Kanhere and R.M, Singru, J, Phys, F. Metal 
Phys. 5, 1146 (19753. 



30 


2) To apply the RFA model to calculate CP and ACPAR 
curves for Ki and Cu for vdilch experimental data as well 
as band structure calcvilations of Et4D exist. An overall 
comparative study of all these results can thus become 
instructive. 

3 ) To test the sensitivity of the REA resxilts to the 

outer electronic conf igxrration (i.e, 3d® 4s^ 3d^ 4s^ 

for Hi) of the transition metals. 

2,2 RFA I4odel ; 

In the RFA model [ 1 ] one starts with the f ree-atom 
Kartree-Fock wavef unctions for outer electrons and truncates 
them at the radius of the Wigner-Seitz (ws) sphere. 

The wavef unctions are then renormal is ed to unity within the 
WS sphere to preserve charge neutrality. In this way the 
free atoms are 'prepared * in an approximate way to enter 
the solid and in a 3d-transition metal such renormalised- 
f ree-atom wavefunctions shoxild give a fair description [ 2] 
of the true ciri^'stal wavefunction at k = O. It has been 
shown [1 ] that the RFA model gives fairly accurate estimates 
of some of the band structure properties for Id-transition 
metals. 

In our calculations we anployed the f ree-atom orbital 
wavef \inctions given by Herman and Skillman [S-J, 

Only the free-4s-abom wavefunctions were renormalised after 
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which they seeme<3 to satisfy the WS boimdary condition 
P^g(r) = O at r = only approximately. It must be 
reraenbered that in this approximation (RFA model) we have 
assigned the momenttm eigenfunctions calcx^lated for the 
k = O electrons to all the 4s electrons. The 3d wavef unctions 
were not renormalised (see Sec, 2,6), 

2,3 Compton Profiles ; 

The procediare follovred by us for computing the CP 
curves was very similar to the one described by Berggren [2], 
For the sake of oampletenessA we reproduce below the in^rtant 
formulas used by us. Ihe contribution by 4s electrons ( k = 0> 
was evaluated by writing the EMD per atom as 

n 

where the sximnation is over all occx:pied states and the 
momentm transform 

A (p) = / d^r e-^P*^ ’J»^(r) ■ (2.2) 

is related to ^ ) by the identity * *^(p) = ) thro-ugh 

o ri ^ T>- o 

^ ^ ^ ^ 
the relation p -- k = 

Before r^lacing the region of integration in Eq, (2,2) 
by the WS inhere an appropriate correction [ 6] was applied. 

The final expression used for the isotropic Oampton profile 
was. 
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J4S (P^l = 4^ I (2,3) 

n=o 

where the f xinctions are given by Berggren [ 2-4 . 

The isotropic contribution by core (ls^2s^2p^3s^3p^3(3.^^) 
electrons was obtained by extrane tight binding approximation 
in which the HMD is written as, 

Pnl^P^ = - jl(pr)P (2.4) 

o 

where the free— atom orbital wavef unctions [ 5] 

and j^(pr) is the spherical Bessel function. The whole 
procedure was tested by calculating the Ctompton profile for 
polycrystalline V, These calculations gave", good agreonent 
with the results of Berggren [ 2 3 for V. 

2 , 4 Positron Wavef Tjnc tip n ; 

The positron vravef unction 4 >, (r) was calculated by a 
plane wave expansion 

I ^ (2.5) 

+ K t 

“K 

described in earlier reports where the resiflts for f (r) 

“T” 

for three symmetry directions [100], [110 J and [ 111] for 
Ni [?] and Cu [8 ] are also shown. For the computation of 
isotropic P(p*),the positron wavef xmction was spherically 
averaged by the method described earlier [?]. 
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2 . 5 Positron Anqijlar Correlation : 

The momenttsn density p(p*) of the tv/o photon pair 
arising from positron annihilation with 4s electrons is 
given [ 6 ] by 

P(p’) = 2 ^ I ^2.6) 

k p’-k, K k=o 
n 


where the sxjmmation includes all the occxjpied states and C . 

k 

the Fourier transform of the product of electron and position 
wavef unctions, is defined as 

Cjp^) = (2tr)''^/2 W ^ (r) (2.7) 

it i + t 


As before, making use of the relation p - Ic = and remembering 

that in the RFA model used by us k = O for the 4s-electrons 

we have (;^ = C (K ) , 
yon 

iC 


Taking a spherical average and integrating as before (2.3) 
vre can get the expression for ACPAR curve for long slit geo- 
metry as. 


N, (p’ ) 
4s z 


= 4w 


I I 


n=o 


C (K )| 
on' 


G (p‘) 
n z 


( 2 . 8 ) 


The isotropic contribution from the core electrons was 

obtained (Compare Eq. (2.4)) by anploying, 

■* ' 




(p*) = (^ I / dr r^ R^(r) P^l (r) j 3 ^{p*r) f 


(2.9) 
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where R_^(r)/r described the spherically averaged positron 
wavef Tjnctions and other quantities vrere the same as that of 
Eq, ( 2 . 4 ). 

2,6 Results and Discussion ; 

2.6.1 Compton Profiles : 

The CP curves for polycrystalline Ni and Cu resulting 
from our calculations are presented in Figs. 2,1 and 2.2 
where they are also compared v;ith the experimental gamma ray 
Compton profiles of Eisehberger and Reed [ 9 '^by normalising 
to ^me total areas (including the extrapolated areas beyond 
p^ = 5.0). It is seen that the RPA model used by us predicts 
satisfactory CP curves for Ni and Cu, as it did for Ti [ 3] , 

V [2 ] and Pe [ 4 1 A band calculation of the CP curves for 
Ki is reported by Rath et al. [ 10] and their curve gives a 
closer agreement with the ejqjeriraent (see Fig. 1, Ref, lO). 

The sensitivity of CP curves (in the low momentum 

region) to the choice of electronic configuration has already . 

been demonstrated for V [2] and Ti [3] , We have r^eated such 

an exercise for nickel, for x-rhich the atomic configuration is 
8 2 

3d 4s , while the configuration compatible with solid state 

9 1 

experience is closer to 3d 4s , Fig. 2,3 shows a comparison of 
the CP curves calculated in different ways. As expected, the 

8 9 

free— atom wavefunctions foif a 3d 4s“ configuration show the 
widest divergence from the experimental curve while the RPA 








CP- Bond profile for Ni 
^ 4S^ free- atom 
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wavef tinctions for a 3d' 4s configuration show closest 
agreement. Earlier APT4 calculations [ 11] for nickel have 
suggested an electronic configuration 3d 4s ^ with 

0,5 < X <1,0, Our calculations indicate that the RFA curve 
based on such a configuration v/oxild give a closer agreement 
with experiment. 

It has been suggested [2] that a better vjay to compare 
the theoretical ciurve with experiment is to subject the theo- 
retical curve to a convolution and deconvolution by the proper 
detector response function. Such a procedure has the effect of 
smoothing the kink in the theoretical curve (Fig. 2,l) and 
improving its agreement with experiment. The RFA curves in 
Figs. 2,1 and 2.2 woxild show a better agreement if such a 
procedure is followed. 

2,6,2 Positron Annihilation : 

The results of our calculations on the momentum density 

for positron annihilation are sinunarised in Tables 2,1 and 2,2 

which show the values of Fourier components C_^ (^<1« (2,7)) 

k 

for the 4s-electrons of nickel and copper. In the case of 
copper a comparison with earlier results (Table 2,2) is 
shovm, the difference between the tw being due to the diffe- 
rent approaches used to describe the electron and positron 


states 
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Table 2.1 


Table of Fourier Components for nickel 


(27t/a) 

[000] 

[111] 

[200] 

[220] 

[311] 

Co (K^) 

(Present work) 

0.925 

-0.078 

-0.075 

-0.034 

+0.008 

Table 

Table 2. 2 

of Fourier Components for < 

Dopper 


(2-K/a) 

[000] 

[111] 

[200] 

[ 220] 

[311] 

Present work • 

0.925 

-0.078 

-0.068 

-0.013 

+0.007 

(Berko and 

0.953 

-0.127 

-0.087 

-0.009 

+0.006 


Plaskett 

[ 6 } 




The theoretical ACi?i.R curves (for long slit geometry) 
for polycry stall ine nickel and copper calculated by us are 
presented in Fig, 2.4 and 2,5 where a comparison with the 
experimental curves and with the band structure calciiLations is 
also given. It is evident that the RPA model does not describe 
the positron annihilation data as well as it does the Compton 
scattering data. In particular the RPA c\irve shows a marked 
deviation (e.g. the deep kink) in the region around the Fermi 
momentum (pp) as well as in the very low momentum region. It 
must be rsmembered, however, that as compared to Compton 
scattering, positrons (because of their positive charge) 
sample the outer electrons with a greater degree of overlap. 

For example, for copper the ratio ^total 

is about 0,5, As a resxilt, the band structure effects assune 
greater importance for positron annihilation. This argtment is 
supported by the relatively better agreement betvjeen the experimoit 
and band structure calculations (shown in Figs. 2, 4 and 2.5), 

This is not stirprising because the band structure calculations 
take into accovtnt the exchange-correlation effects of the 
d-electrons, hybridization and most important of them sill the 
details of the Fermi sxorface shapes. The residual disagreement 
(Pigs. 2,4 and 2,5) between the experiment (curve b) and the band 
structure calcTjlations (cxirve C) coxold be due to possible 
enhancement effects, uncertainties in the spherical averaging etc. 
However outr results suggest that RFA np_dcl gives a 'poof . 
description of the 4s-eloGtrDns vis- a- vis positron annihilation. 
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It might be argued that in viev; of the large nixnber 
of 3d electrons in nickel and copper the effect of renormali- 
zation of the 3d— states should be incl’jded [2], Our calctJla- 
tions indicated that renormalization of Id-states had the 
effect of broadening the curves and v/idening the deviations 
between RFA theory and experiment. It may be pointed out that 
a comparison of different results of positron annihilation 
studies in copper [12 ] indicates that the Id-electrons are 
more spreadout in position space than the free-atom wave- 
functions indicate. This might perhaps be due to the possibi- 
lity that the effect of renormalization of 3 d-wavef unctions and 
the effect of the d— hybridization with conduction band act in 
opposite directions. Considering all this, we did not renormalise 
the Id-states in our calculations, 

2.7 Conclusions : 

Referring bade to the motivation behind this work we 
conclude that 

a) RFA model cannot give a satisfactory description of 
the positron annihilation process for which a detailed 
calculation anploying band structure is necessary. 

b) It gives a satisfactory but gross description of the 
overall features of the Compton profile of Ni and Cu, 
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cO Our results sij^sport an electronic configuration closer 

to 3d^ 4s^ for Ni. 

However, the finer details and particularly the 
anisotropies of the Compton profiles canrot be described 
without a proper band structure calculation and hence the 
applicaioility of the RFA model is limited. 

In the next Chapter v/e shall describe our band structiure 


calculations of EMD, 
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Bmo simja-m mrm im calcuiation of 
mmmm pensities 

3,1 Introduction : 

It was seen in the last chapter that for a detailed 
and quantitative analysis of the CP and ACPAR data^ especially 
for predicting the anisotropies, theoretical calcxslations 
aPEiploying band structure methods are essential. For doing 
band structure calculations one has a choice of using any 
of the available standard methods (e.g. APW, KKR etc.) vhich 
are known to give good results for other ground state 
properties of transition metals. As emphasized in Chapter 1, 
our aim is to analyze the EMD in a physically transparent 
manner, clearly bringing out the role of d bands, the 
structure of Fermi surface and in general the electronic 
structure of transition metals, Iherefore, instead of opting 
for an accxirate, sophisticated and time-consuning method, we 
tried to seek a fast, but approximate method which can meet 
oixr needs. lb wards this aim we have chosen Hxibbard's 
method [l] especially developed for transition metals. Ihis 
method has already been employed to calculate the IMD in 

V ■ 

bu [2,3] Fe [2,4] and Ki [5], 

In the following we first give a brief outline of the 
Hubbard's method and obtain an expression for ■'■tc aiD, 

p (^, following Mijnarends [ 2 ] . More details of the method 



47 


can be foimd in literature [1/2], The results obtained by us 
by this method for Pd and Ag vd.ll be presented in the next 
chapter, 

3,2 The ^ubbard-Miinarends Method ; 


In the KKRZ formulation, Ziman’s determinantal equation 

» « 

giving the energy levels £ for the electron in state k can be 
v/ritten 1 as follows: 

(l4 _c) + E r E^, = 0 (3.1) 

n* nn 

In 2 »bove = k + K^, the being reciprocal lattice vectors? 

are the coefficients of pseudo wavef unctions. The pseudo- 
potential matrix elements are given by. 


^nn* ^ 1 ^1 ^1 ^1^ P^^Ccos € 

(3.2) 


nn’ 


where. 


tan n, 

^ = - : 

(3.3) 

Here f is the volime of atomic polyhedron, E* = ^^(E) 

is the phase shift of the 1-th partial wave scattered from 
atomic polyhedron (mxiffin— tin potential) at energy E, and 
n., are the spherical Bessel and Nermann functions respectively. 
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e . is the angle between k and k ,, P, is the Legendre 
nn n n 1 

polynomial of order 1. One can choose r^ arbitrarily in 
the interval (O, being the radius of the muff in- tin 

sphere. 

In the case of the simple metals, tann^(E) is a 
smoothly varying function and is small. However, in the 
case of transition metals (i.e. d bands), tan 1^2 ajid hence 
also Y 2 passes through a narrow resonance at low energy. 

In general v/ill be small, except near the resonance 

in tan One can then write 

^1 ~ ^1 “ (3. 4 ) 


where 7 is the width of the resonance occurring for the 
partial wave 1 =!XatE = Now is always snail, 

since the resonant part is siibtracted out. We can rewrite (3,l) 
in the following form [l]. 






(3.5) 


where the matrix elements of the siibrnatrices are given by, 


K 


hn* 


n nn' 


{3,6a) 


V 


nn' 


~ r (21+1) r^^Yj ji(\ ®nn*^ 


(3. 6b) 
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= -4fr. j (Ic r ■) Y, (5 ) 

rm X T -^X n^Am'n 


(3« 6c5 




EQ-e T 


(3.6d) 


The matrix |^ is of the order (2^ + l); matrix Q contains 
the non-resonant part and the resonant part of the 
wavef unction. 

The above equation (3.5) can be reduced to one with 
a lower order by separating the vectors in two sets, a 
preferred set P and the r^aainder R. Introducing the approxi- 
mations \//p^ = \^p = ^ andxGEneglecting with respec .? 

to (!^ - e ) when belongs to the set R, one can get, 

V h„ A-«i/la; 


- e 


L + h„a. = o 

R R 


(3.8) 


where 


, _ - (4'g) ^..2 (r) j ix^^ifx?. ^ -iv* ft \ r 

rom ' mm T^Xr^^ ^_2 Xm X ro * '^n ^ 


^^-n - ^ 




(3.9) 


(r) 

and where the sot S is over all the vectors of the 



50 


* remainder’ set R, For the case of d-reso nance the order of 
Eq, ( 3 , 7 ) is (n^ +5) vrhere is the nimber of preferred 
vectors. 

The form of (3,7} is very similar to the one used in 
model psexKiopotential schemes [6 ] , The suianatrix describes 
the nearly free electron (NFE) part of the conduction band (S-p) 
while the 5x5 submatrix describes the d bands and the 
off-diagonal matrices and ^ cause the hybridization 
between these bands. The sufcsnatrices ^ and are energy 
dependent and the next step is to remove the energy dependence 
in the matrix n so that the problem is reduced to that of a 
di agonal ization. This can be achieved by suitably approximating 
the energy dependence of ^ and V as follows# 


^ / 5< + y/ck^) b 




(3.10) 


where is the mean d band energy. 

A nimber of refinements of above approximations are ■ 
suggested by Hiibbard [1 ] . In essence the problem rediKies to the 
solution: of the secular equation. 




o 


(3.11) 
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where gives the rson-resonant part of the wavefunction 

while gives the resonant part. Ihe remaining coefficients 

can be determined v/ith the help of Eq. (3,8), 

I>loyd J has shown that in the region outside the 
muff in- tin sphere (radius r^) the wavefunction can be ejqsanded 
in terms of plane v/aves having as the plane wave coefficients 

^ _ z 

h (t)1/2 ^ 

For r < i.e. within the inscribed sphere the wavefunction 
can be ejqianded in terns of spherical harmonics, 

. A mr n F-I ( r ) ->■ 

(r) = — J h — = r (rf (3.13) 

k (T)1/2 L ^ ^ 

where the coefficients Xj^ (where L = 1, m) are determined 
by the contin-ulty condition at r = r^, i.e. 

3l'-n <V 

n 

In practice only the first fevr 1-values are retained 
in (3.13) and (3.14) and this truncation results in a slight 
discontinuity at r = r^ , 

In the above we have sketched a bare outline of 
Hubbard’s theoiry details of v/hich can be found in the 
literature [!]• fhe next step now is to extend this method 
for calculation of electronic wavefunctions and momentum 

C£Ni H , -.i.- 

hoc* No. i''h. ® 


B 0 


i (^+^) .r" 


(3,12) 


density 
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In the independent particle model the electron momentian 
density is given by^ 

P (p) = I f(3t^j) {a. (ic/ 'p)5^ 6(p - it - K . ) (3,15) 

where f(k, j) is the Fermi-Dirac distribution function such that 
at T = O K, f (k, j) = O if the state (k^ j) is unoccupied and 
f (k, j) =1 if the state is occupied. In (3.15) j is the band 
index while denotes the reciprocal lattice vector, where 

(k, p) is the Fourier" transform of the electron wavefunction, 
and is given by, 

A. dC p) = / e“^P’^ ^ (r) dr (3.16) 

^ crystal k, j 


The above expression can be simplified by using Eq. (3.12) 
and (3.13) for the electron wavef unction. Using these 
equations, Eq, (3,16) can be written as. 


Aj (k, p) 





jldVPi 

l’'i. - p1 


^ (4t)^ 


3 

S E 
1=0 m R^(r^) 



r . 

-1 •? 

X Y^jj^(p) / j^(pr) R^(r) r dr 
o 


(3.17) 


Equations (3,15) arsi (3.17) form the basis of our computation 
of momentum distribution due to 4d (or 3d) and condtsction 
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electrons. Contribution to the EMD from core electrons 
could be obtainei by employing Eq, (2,4) and is already 
tabulated in the literature [8]. 

3.3 Computational Detail s i 

The computational procedure followed by us was divided 
into following stages (Pig. 3.1). 

a) Calculation of Phase Shifts and Radial Wavef -unctions t 

The starting point for the calculation of the phase 
shift is the muffin-tin potential [9.1 interpolated on an 
appropriate Herman- Skillm an mesh I 10 1, linploying the 
muffin-tin potential, phase shifts were calculated over 
entire energy range of interest at about 200 energy values 
with an interval of 0,01 Ry, (atomic units are used for other 
quantities). The radial part of the wavef unction was* stored 
on a magnetic... tape for further \ase. This stage of the calcu- 
lation was independent of the crystal structure, 

b) CalcTxlation of Energy Eigenvalues Wavef unctions 

(B^ and coefficients) and Momentm Density P(p): 

Employing the phase shifts and radial wavefiinction 
so obtained, we next calculated E^, and P(p) 

(see Eq. (3,l) to (3.17)) for each selected k-point. 

Partial waves r^jto 1 =3 were inclx^ded in the expansion 

(3.13) of the wavef unction inside the inscribed inhere while 
89 reciprocal lattice vectors were used in the e>qpansion 
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(3,12) outside the sphere, Mijnarends [2] has reported a 3% 
increase in the p for Fe in the intermediate nKjmentim 
region due to inclusion of 1=3 term. We have found that 
inclusion of 113 plane waves (instead of 89} had marginal 
effect on P , 

The accixracy of Hubbard's method has been discussed 
by Kaga [ll ] for the case of Pd and he has pointed out that 
the major errors in the energy values arise from the linear 
inteirpolations and extrapolations of the matrices /i\ (e) and 
(e) and from the small nrsnber of plane waves in the ‘pre- 
ferred set* P, We have studied [12] the effect of increasing 
the number of plane waves (in the P set) from 4 to 14 for Pd 
and Ag, Such an effect on energy values and momentum density 
P(^) is shown in Tables 3,1 arsi 3,2 at some syirmetry points. 

It can be seen that the energy levels of Pd near the Fermi 

energy (E ) are fairly stable and are not affected significantly 

i? 

but energy values of the lower bands increase by as much as 
0,03 Ry, bringing them in closer agreemoit with other results [13] 
as shown in Table 3,3, The effect of increasing the *P' vectors 
is greater on P(^, thus> for example at (p = O^ P(p) 
decreases by as much as 20% while the energy eigenvalue at 
increases by 5% only. In the case of Ag ( .Table 3.4 ) 
we found that the maximum change in energy was of the order 
of 0,01 Ry, while the change in P(p) at ^ = O was only about 
1%. In view of these observations we employe! only 4 vectors 
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Table 3,1 

|5 

u 

Effect of the n\snber of preferred vectors on the energyvalues 
of Pd, The first colimn gives k vectors in units of 2i5/a 
alongwith the BSX'J symbol. The second column specifies the 
nxmber of preferred vectors, and the remaining colxmns give 
the energy values in Rydberg, for the lowest six bands. 


1- 

No. of 


^ENERGY IN 

Ry. 





i 

2 

3 

4 

5 

6 

r 

(0, O, 0) 

4 

0.104 

0.304 

0.306 

0.306 

0.425 

O. 427 


14 

0.101 

0.319 

0.320 

0.320 

0.427 

0,427 

X 

o 

H 

•• 

o 

4 . 

0.107 

0.117 

0.5C7 

0.536 

0.536 

1-.141 


14 

0.131 

0.161 

0.505 

0.533 

0,536 

O. 855 

K 

(0.75,0.75,0) 

4 

0.156 

0,218 

0.387 

0,447 

0.505 

1.236 


14 

0.179 

0.220 

0.339 

0,446 

0.503 

1.105 

L 

(0.5, 0.5, 0,5) 

4 

0.139 

O. 296 

0.29G 

0.511 

0.511 

1.194 


14 

0.159 

0.301 

0.312 

0,510 

0.511 

O, 652 

n 

(0.25, 1.0, 0.25) 4 

0.157 

0.215 

0,339 

0.447 

0.508 

1.111 


14 

0.130 

0.227 

0.390 

0. 443 

0.504 

1.106 

* All values are with respect 

to the 

zero of 

the muffin tin 


potential which is at -0,85007 Ry, 
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Table 3,2 

Effect of the ni^nber of preferred vectors on the momentun 
density p(p) (for p = for Pd, np denotes the nxinber 
of preferred vectors used. 


• 4 * 

n N 

P- 

r 

X 

K 

L 

U 

4 

0.912 

0.197 

0.255‘ 

0.193 

0.280 

14 

0,865 

0.175^ 

0.248 

0.17 4 

0.267 





Table 


Effect of the ntimber of preferr 
of Ag. All other quantities ir 


that for Table 3.1, 
the values of P (p) 


The last c 


:or the two 


No, of 

•p’ 

vectors 1 


r 

(o, O, O) 4 -0.025 

14 -0.025 

X 

(o, 1 , O) 4 -0.067 

14 -0.057 

K 

(0.75,0.75, O) 4 -0.038 

14 -0.029 


4 -0.042 

14 -0.037 


4 0. 023 

14 0.024 
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ENERGY m Ry . 



‘^t 

2 

3 

4 

5 

6 


0.042 

0.04S 

0.043 

0. 04S 

0.043 
O. 048 

0.111 

0.111 

0.112 

0,112 

0* 811 
O. 32- 

•0.065 
■0. 047 

0.160 

0.159 

0.176 

0.176 

0. 176 
0. 176 

0.645 
0. 645 

0.141 

0.131 

0.000 

0.000 

0.101 

0.101 

0.127 

0.126 

0.159 

0.157 

0,818 
0. 817 

0.171 

0.17! 

0.039 

0.039 

0.162 

0.162 

0.477 

o.ssi 

o. 045 

0, 045 

0.162 

0.162 

0.477 

0.5 si 

0.037 

O. 040 

0.085 

0.081 

0.089 

0.086 

0. 089 
0. 089 

0.262 

0.260 

O. 79 I 

0.781 
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Table 3.5 

Some parameters used in the calculations. 

P d [14] Ag [15 ] 

a, lattice constant (a, u. ) 7.3530 7.7112 

Constant potential between the 

spheres^ to be taken as Muffin- -0.85007 -1,1183 

Tin zero (Ry.) 

Ej,, Fermi energy with respect 

to the zero of the Muff in- Tin 0,5246 0.518 

potential (Ry.) 

rj^/ Muff in- Tin radius (a.u, ) 2.59968 2.7 26 

e the mean d band energy 

(the energy at which tan ^2 0.335 0,070 

goes through a resonance) 

* The values quoted are from our calcifLations. 



61 


Table 3,6 

The preferred vectors (in units of 2Tr/a) used in the calcula- 
tions for Pd (14). In the case of Ag only first fo'or vectors 

are used* 


1 

(0, 

0, 

0) 

2 

(-1, 

-1 

. -1 

3 

(0. 

-2. 

0) 

4 

(-1, 

-1 

. 1) 

5 

(1. 

1. 

1) 

6 

(1/ 

t,- 

-1) 

7 

Cl. 

-1. 

1) 

S 

(1. 

-1. 

-1) 

9 

c-i. 


1) 

lO 

(-1. 

' 1 % 

-1) 

11 

(2, 

0, 

o) 

12 

(-2, 

0, 

0) 

13 

(0, 

2/ 

0) 

14 

(0, 

0. 

2) 
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in the ’P' set for Ag (thus saving computational time) 
while 14 vectors were used for Pd. Other numerical data about 
the parameters employed in our calculations are given in 
Table 3,5, v;hile Table 3.6 lists the set of preferred vectors 
used by us. Alxjut 90 k-points in 1/48 of the Brillouin Zone 
(in the (lOO) and (lio) planes) were used to calculate the 
energy eigenvalues. Employing suitable K-vectors, P(p) was 
calculated for about 800 p— points in the above symmetry planes. 
Typical computer time taken to calculate .^ID at one k-point 
(along with its U— processes) vras about 2.0 minutes on the 
IBM 7044 Computer of the Indian Institute of Technology, Kanpur. 

The results of our calculations for Pd and Ag are 
discussed in the next chapter. 
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ChapteA. 4 


ELECTROM fiOACMTUW VEf^SITV T« Pd* AM? % 


4. 1 Introduction : 

In this chapter we present oxir results for the 
band structure and IMD for Pd and Ag and discuss their 
important features. Our restilts for the energy bands and 
Fermi surface are presented and compared with other theore- 
tical results in Sec. 4.2 v/hile the results for total UdD 
as well as the band-by— band contributions to the total IMD 
are discussed in Sec. 4,3, In the Sec. 4.4 the HMD of Pd 
and Ag are compared with that of Ni and Cu. 

4. 2 Energy Bands and Fermi Surface of Pd ard. Ag : 

Out resiilts for the band structure for Pd and Ag 
along important symmetiri'' directions are shown in "‘Fig. 4.1 
and 4,2 respectively. For Pd, we have employed the muff in— tin 
potential which was constructed by Wong [1] using a non- 
relativistic free-atom HFS wavef unctions for a 4d^^5s*^ 
configuration. The agreement between our resiilts (Fig. 4,l) 

* Preliminary results of our calculations for Pd are 
published in the paper: D.G. Kanhere and R;M. Singru, 

Phys. Letts. 53A , 67 (1975), 
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and those of Wong [ 1] is good, A detailed sttKJy of the 
electronic structure of Pd is made by Mueller et al, [ 2], 

These authors have employed an APW method, as well as a 
combined interpolation scheme, and have included relativistic 
effects also. Kxogh- Anderson [ 3] and Tex-mri [4] have also 
calculated the band structure of Pd, Our calculated 
eigenvalues for some ■ symmetry points for Pd are compared 
with those of Mueller et.al. [2] in Table 4,1. In general 
the agreement of our results with other restalts [1,2 ] is 
good. The width of the d band (X^ — X^') observed by us is 
closer to the results of Mueller et al. [ 2 ], The cross 
sections of the Fermi surface for Pd by the (lOO} and (llO) 
plane are shown in Pig. 4.3, The Fermi surface of Pd 
consists of (i) a T -centered sixth band electron sheet 
which displays large distortions. This sheet is primarily 
d-like, (ii) X»c entered hole sheet formed from the fifth-band 
electrons giving rise to an open netvjork in the extended zone, 
(iii) snail hole pockets around X, Above featxnres of the 
Fermi surface (Fig. 4.3) agree qualitatively with other 
calculations [2—4] and experimental results [5, 6|, A L— centered 
fifth band hole pocket found by Mueller et al. [ 2 ] is not yet 
observed ejqjerimen tally and is absent in our resiilts. This 
difference could perhaps come from the relativistic effects. 

Band structure of Ag has been calculated by Bhatnagar 
employing APW method, by Snow [s ] anploying self-consistent 
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Table 4,1 

Comparison of some important band widths for Pd, obtained 
by us with others. All the energy values are in Rydbergs. 


Present Mueller et.al 
worlc [ 2] 


VJidth of the (sp) band 

Width of the H’band 

Separation 

Distance of Ep to 
upper most *d' level 


i 


^4 

- "i 

0,754 

^5 

- ^ 

0,404 

< 


0.310 



0,017 

5 

F 



0.751 

0.396 

0.30S 


0.017 


0.018 
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Table 4.2 

Comparison of some iirportant bandwidths for Ag, 
Energy values are in Rydbergs. 


Present Bhatna- Rijsenbrij and 
worlc gar [7 ] Pondse [9j 


Width of the (sp) band 

X - r 

4 1 

0.670 

0.652 

0.688 

Width of the ‘d* band 

X5 - X^ 

0.243 

0.240 

0.251 

S^aration 

^4 - ^5 

0 . 469 

0.457 

0,442 

Distance of E^, to 
upper most 'd* level 

Ep - Xg 

0.33 2' 

0.322 

0,292 
















/ 

APW method, by Rijsenbrij [9] anploying non-relativistic 
KI<R method and by Christensen [lO] by relativistic APW method. 
We onployed the muffin-tin potential for Ag constructed by 
Bhatnagar [ 7 ] . A comparison of our results with other 
calculations is shown in Table 4.2 and it shows an oversill 
satisfactory agreement. Fermi surface of Ag (Fig, 4,4) 

shows features common to other noble metals (viz, Cu and Au) 

♦ 

in that it is in contact with the zone face at L giving rise 
to a neck along eight [3.11] directions [11 ]. 

4,3 Electron Momentun Distribution for Pd and Ag : 

The results of our EMD calculations for 4d and 
conduction electrons in Pd and Ag are presented in Pig. 4. 5 
to 4.24, The EMD for Pd along the important symmetry directions 
A (lOO), r (110) and A (111) is shovm in Pigs, 4.5 — 4,7. 

These curves show the total HID as well as the band— by-band 
EMD with the relevant band structure shown in the inset. 

The labels a,b, correspond to the particular band which 
contributes bo the EKiD while the dashed part of the curve 
indicates the contribution from the unoccupied bands. The 
results for Ag are sho'^'m in Pigs, 4,8 — 4,10 in a similar 
manner. 

The curves describing total EMD show significant 
strucfure (i.e. valleys, steps etc.) which can be attributed 
to the sixth band electron Fermi surface and its reflections 
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X 



in atomic units 




DiO 
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in higher zones (through Unklapp processes). In particular 
prominent steps are observed for Pd along A and A directions. 
A smaller step is seen along the ^ direction and it can be 
explained as follows. The hybridized ^ band (marked as ’c* 
in Pig, 4,6) crosses the Fermi level at a smaller k-value and 
its" v/avefunction retains the original d character at the 
Fermi level. The contribution of this band goes as p"^ and 
hence has smaller value at the Fermi level. The EJID in Pd 
shows strong momentum dependence and this can be ascribed to 
the s-d hybridization. If we examine the contribution to 
the SMD by individual bands# we notice that only the 
and bands alone contribute. This is a resiilt of a 
selection rule arising from a full symmetry possessed by 
P (p) , as pointed out by Mijnarends [12], 

The effect of ;sp--d hybridization on S«!D is shown 
by the individual bands but the net effect of hybridization 
on the total EMD is small. The contribution of the lowest 
A^# and bands (marked *a* in Figs. 4, 5-4,7) falls 
off as k increases because the bands hybridize with the 
higher bands. The contribution due to the higher A and 
bands (marked *b* in Figs. 4.5-4.7) varies initially as 
p^ (as expected for a d band) but rises faster due to the 
mixing of s— character. Along the Z direction# there are two 
higher bands (marked ‘b* and *g‘ in Pig, 4.6) axA among 

these the *b * band is d-like throughout the zone while 'c * 
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hybridizes vrith *a* strongly and their respective EMD slx)W 
this character. A prominent featizre of all the EMD curves 
is the significant value of p (p) at high momenta 
v/hich is contributed by the d— lilce states via the ISnlclapp 
processes, I-fowever, the momentun dependence in the region 
p = 1,0 - 2,5 a,u. along 2 direction is significantly 
different from that along and A directions. This behaviour 
of Er4D is expected to manifest itself in the anisotropy of 
directional CP, 

The behaviour of the EMD in Ag along the three 
S3mimetry directions (Figs. 4.8 - 4,10} can be understood in 
a manner similar to Pd. The overall features of the E14D 
in Ag shov? similarity with that of Cu [12] as expected. 

As Compared to Pd, Ag has narrower d bands which lie v/ell 
belov7 the Fermi level, A.s a resxiLt the states in Ag near 
Fermi level are almost completely hybridized and possess a 
dominant p character. These features of the bands are 
reflected in the sharper (as compared to Pd) steps, 
especially along the 2 direction in the first zone and 
relatively smaller steps in higher zones. The nearly-free 
electron (KFE) character of the in Ag is shown by its 
flat behavio^l^ in the first zone. Remanbering that the 4d 
band is completely full in Ag, this resrdt is not unexpected. 
The effect of the neclc in the Fermi surface of Ag manifests 
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itself in the smooth behaviour of the EMD along the A 
direction, although the slope of EMD at the zone face is 
considerably large. The half-v;idth of the IMD contributed 
by the lov/est bands (marked *a*) are analler in Ag than in 
Pd. This difference may be due to the smaller degree of 
hybridization in Ag, The EMD cur^/es for Ag sliow less aniso- 
tropy in the region p < 1 a.u. when compared to Pd. 

A deeper insight in the EKD is obtained if it is 
plotted as a contour-wise diagram in a plane. Such plots 
are especially useful to study the effects due to (i) anall 
changes in band structiire near the Fermi level or (ii) changes 
in the Fermi level itself due to alloying, which can be 
analyzed in the framework of rigid band model. We have, 
therefore, plotted the isodenses for Pd and Ag in the planes 
(lOO) and (llO) in Figs, 4.11 - 4,24. The plots of total 
EMD contours are followed by the plots of EMD due to invidi— 
dual bands 1,2,..., 6, the numbering of the bands being 
done according to increasing energy (Pig. 4.1 and 4,2). The 
dotted curves (Figs. 4.11 to 4.24) represent the lines along 
which the bands touch. A transfer of EMD takes place across 
these lines as the symmetric {+ or -) character of the band 
changes. As a result of the selection rule [l2] discussed 
above only the bands with (+) symmetry can contribute to the 
EMD. In Figs. 4.15, 4.16, 4.17 and 4.24 we have also shown 
the Fermi surface by indicating the vroccrpied region by 
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shaded The +/- signs indicate the representation of 

the hand electron wavef ’unction in that region. The zone 
iDoundaries are shovm by thin lines whereas the isodenses 
(con-cours) are shov.m by thick curves and magnitude of 
the a-'lD is indicated by the nimbers 0.05, O. i, 0,2, 0,8, 

shown on the contours. 

The total EllD for Pd (Pig. 4.1l) in the (lOO) and 
(11 O') planes shows strong anisotropies even in the first 
zone. The isodenses for P (^ = 0,3 - 0,5 show bxilges in 
the (lOO) and (liO) planes, arising mainly from the first 
band and partly from the second band (Figs. 4,12 and 4,13). 
This result, indicates that for Pd even the low-lying, 
completely filled bands contribute to the anisotropy of the 
Q^ID, This aspect is made more clear if the isodenses of 
band 1,2 and 3 (Figs. 4.12 to 4.14) are closely examined. A 
study of the isodenses in higher zones indicates that their 
contribution comes mainly from the d bands via the Unklapp 
process. 

In the (llO) plane, band 1 and 2 contribute all over 
the first zone because their v/ave functions have (+) symmetry 
and the selection rule [12] allov/s their contributions. On 
the other hand, the remaining bands contribute only in 
those regions of the zones where they have proper (+) 
symmetry. The validity of this selection rule in our 
results is nicely brought out by the isodenses. The remarks 



made earlier about the effect of S/P/d... character of the 
bands on the EMD can also explain the nature of &iD in the 
(lOO) and (lio) planes. Looking at the isodenses of band 4 
(Fig. 4.15) we notice that the hole pockets centered arotind 
X cannot be observed by Compton profile e>g3eriments because 
of the negligible amount of EMD surrounding them, Ihe fifth- 
band X— centered hole sheet (Jungle-gym structure) is unobser- 
vable in the (llO) plane owing to the wrong (-) symmetry of 
the fiftli band. Examination of Pig. 4.17 st^ows that the 
sixth band contributes to the EMD mostly in the unoccupied 
parts of the zone. The effect of the T— centered sixth 
band Femi surface on EMD is most significant along A and 
E directions but it becomes less significant as one moves 
away from these directions. This behaviotir is expected to 
manifest itself in the anisotropy of Compton profiles. It 
is interesting to note that for the sixth band the contribu- 
tions in higher occupied zones due to the n-processes are of 
the same order of magnitude as that in the first occupied 
zones. 

The isodenses for Ag can be analyzed in a manner 
similar to the discussions given for Pd. Generally the. 
anisotropy observed for* Ag is less than that observed for 
Pd. The total EMD (Fig. 4.18) in the (lOO) plane falls off 
rather sharply in the first zone but the decline is slower 
in the higher zones. This behaviour is larg^y due to 
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the fact that fox* Ag the Bth band, contributes significantly 
to the E2-iD near Fermi surface. The effect of the necTc of the 
Fermi surface of Ag is clearly visible along the [ill] 
direction near L in Figs. 4.18 and 4.24. For F^g^ band 1 
contributes throughout the first zone in the (lio) plane, 
but in the (lOO) plane it contributes only in those parts 
where it has (+) character. The contribution of the sixth 
band in the occupied part of the higher zone is much smaller 
than that in the first zone, A comparison of the total EMD 
contours (Fig. 4.11 and 4.18) shows that in the first zone, 
the effects of the sixth band Fermi surface is more important 
for Ag than for Pd, 

We feel that our results for the total EMD (Fig. 4.11 
and 4.13) for Pd and Ag can be useful for analyzing a future 
unfolding of experimental CP data for the reconstruction of 
St4D. The isodenses for individual bands, however, canngt 
provide a direct comparison with experimental CP data which 
are built only by the total EMD. Nevertheless, they can be 
compared with other theoretical band structure calculations 
of EMD to examine the effects of band structure on IMD, It 
may be pointed out that, in principle, the method of quasi- 
free (e, 2e) scattering studies can measiare [13, 14] the 
energy depencient momentum densities E)* 
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4.4 C omparis on o f EM D: Ni Pd^ Cu vs Ag. 

A comparison of the total EMD for the isoelectronic 
elements Ni [15] and Pd along the symmetry directions [lOO], 

[ lio] and [lll]shov7s that although the overall features of 
the EMD are similar, there are significant differences in 
their EMD along [ill ] direction. These differences are 
due to the 6th band Fermi surface. It must be remembered 
that Ni being ferromagnetic a direct comparison with Pd is 
not possible. tibwever, the minor differences along [lOO] 
and [llO] direction can be understood in terms of the 
crossing of the sixth band with the Fermi level. 

A comparison between the isoelectronic elements 
Cu [ 12 ] and Ag reveals a ranarlcable similarity in their EMD, 
an observation to be esqpected from their similar band structure 
(characteristic of noble metals). Our results show that 
small changes in the' neck radii (along the [111] direction) 
of these t'wo metals do not show up significantly in their 
total EMD. 

It would be interesting to compare the results of 
EMD for ACFARand CP employing same electron wavef lanctions. 

Such results for Pd [ 16 ] indicate a similarity between the 
two EMD' s in the first zone. However, differences appear 
in the higher zones as CP shows significant high momentum 
components. As pointed out in Chapter 1, this behaviour is 
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to be expected because the positron is bept away from the 
nucleus and annihilates preferentially v/ith the outer 
(conduction) electrons. 

4 , 5 Lixai tatlons p_f _ Our Calculat ions: 

In the end we wish to point out some limitations 
of our calculations, 

i) Firstly the approximations inherent in the Hubbard 
scheme are expected to limit the accuracy of our wavef unctions. 

ii) Our calcuiations do not take into account the relati- 
vistic effects which can be inportant for 4d metals. Although 
the qualitative features of the EMD for Pd and Ag are not 
expected to be affected much by relativistic effects/ a 
quantitative estimate is not yet available, 

iii) We have calculated the EMD only in the two planes at 

limited number (“ 90) of k-points. We have not calculated 
the EMD beyond p > 3.5 a.u. because of the limitations of 
our metliod. This would limit the use of our results for the 
calculation of the directional profiles using (1*9 ). 

Inspite of these limitations we feel that our 
results for Pd and Ag have given a good understanding of 
the EMD in these metal S/ underlining the importance of band 
strucfure methods for transition metals. The results of 
the present chapter become more directly applicable when 
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the methods outlined in Chapter 6 are employed to calcxilate 
directional profiles and, anisotropies. The resiilts suggest 
that experimental study of CP in Pd and Ag single crystals 
should be undertal^en and our res’olts (Chapter 4 and 6) 
can then provide a basis for a proper comparison between 
theory and experiment. It is known that Pd forms interest- 
ing alloys v^ith Ag involving interesting changes in the 
electronic structure. Such Pd— Ag alloys are being studied 
by ACPAR [ 17] . We feel that o\ir results for the EMD of pure 
Pd and Ag can provide an introductory basis to understand 
such studies and in particular to examine the role of 
positron in ACBiR, 
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ChapteA 5 


S!/STEMATIC$ OF THE ELECTWM MOMENTUM 
VEUSITV IW TRAtiSinOM METALS 

5,1 Introduction ; 

Because of the interesting electronic structxire of 
the transition metals, their energy band structure and 
other ground state properties have been a subject of great 
theoretical interest. Snow and Waber [ 1 ] and Matheiss [ 2 ] 
have systematically studied the band structure and other 
properties of the 3d metals, Ti through Cu, Although some 
band structure calcifLations of the EMD (relevant to CP as 
well as ACPAr) have been reported in literature, each of 
these studies pertain either to one or two metals. Similarly 
different band structure methods have been used for different 
metals. A systematic study of the EMD for the iron-series 
3d metals employing the same band structure method can 
therefore be interesting. In this chapter we present our 
results for the EMD (pertaining to CP) for the 3d metals 
V, Cr, Pe, Ni and Cu calculated by a:nploying Hubbard-MIijnarends. 
method. These results are then qualitatively discussed in 
terms of their respective band structure and systematic trends 
are pointed out. 
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5*2 The Band. Structure of 3d Metals : 

The band structure and the EMD for V, Cr, Fe, Ni and 
Cu were calculated along the three symmetry directions A t 
2 ; and A by the method described in Chapter 3, The muffin— tin 
potential for Cr was generated by Lovrdin*s cc —expansion 
method [ 1/3] using full slater exchange. Atomic charge 
densities for the configuration 3d^ were generated by 
using Harm an- Ski 11m an programme [4 ], Muffin— tin potentials 
as reported by Mijnarends [5] , Wood [6]/ Walcoh [7/8] and 
Biardick [ 9 ] were used for V/ Fe, Ni and Cu respectively. 

Other computational details are common and they are already 
discussed in Chapter 3, 

TaJDle 5,1 shows the lattice constant/ important 
energy bandwidths/ location of the resonance S ^ for the d 
phase shifts and location of Fermi energy (% ) with respect 
to the bottom of the conduction band ( r^) for these metals, 

A detailed discussion of the systanatics displayed by the 
energy bands etc. by these metals is given by Snov7 and 
Waber [l]. Our present calculations show trends similar 
to these. The band structure of all the transition metals 
is characterised by the presence of a fairly narrow d band 
that overlaps and hybridizes with a broader (and nearly 
free— electron) s— p band. As one goes along the iron 3d 
series one observes the following characteristics? 
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Table 5,1 


Some parameters and bandwidths for the 3d transition metals. 
All the energy values are in Rydbergs. 



V 

Cr 

Fe 

■Ni 

Cu 

'I 

a/ lattice constant (a.u. ) 

5.7225 

5,4494 

5.406 

6* 65 

6.8309 

Cd 

0.868 

0,892 

0.492 

0,534 

0.306 

^ 

0.54 

0.58 

0.59 

0,67 

0.66 

Width of the 'd' band 

0.55 

0.52 

0.40 

0.36 

0.265 

Width of the occupied 






‘ d ‘ band 

0,23 

0.25 

0.35 

0.36 

Compl etely 
occupied 


where 




^2 

Ep - X, = 
®P - ^12= 


Location of the resonance in tany'j ^2 

Position of the Fermi level with respect to 
the bottom of the conduction band. 

Width of the 'd' band for fee metals 

Width of the 'd’ band for bcc metals 

VJidth of the occupied ’d' band for fee metals 

Width of the occupied ’d' band for bcc metals. 
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a) As the atomic number increases (Z = 23 to 29) the 3d 

states are grouped into small energy range, i.e. the d-lilce 
states are more and more tightly bound, Ihe width of the 
total d band (X^ - or shown in Table 5,l) 

decreases from V to Cu, but the width of the occupied d 
band or E^, - increases. 

b) Compared with the bottom of the conduction band, 

(Ep •* Permi level (Ep) rises with the atomic number 

because more and more d electrons are being accommodated, 

5.3 ResxfLts for the EMD of 3d Metals and Discussion ? 

The plots of the total as V7ell as band-by-band EMD 
for V, Cr, Fe (all bcc), and Ni and Cu (both fee) along the 
three symmetry direction a / Z and A are shown in 
Figs. 5.1 - 5,6, The relevant band structure is also shown 
in the figures. It is more convenient to divide the discus- 
sion about the systematics of these resiilts according to the 
symmetry directions. 

a) A ( 100 ) Direction t 

The striking feature displayed by the total EMD is 
the sharp breaks occurring for Fe, Ni and Cu at their 
Fermi momenta (pp) . In the case of V and Cr only the 
occupied lower A^ band contributes and the EMD shows a 
smooth character with no breaks at Pp. ^ small step appears 
for Pe as the higher band becomes partially occx:^ied. 





Pig* 5i*2« Bn«r«y bwoida and f(p) for and Cu along 

riOOJ^direotion# ©it ▼trtleal aoala of taexgy 
intthd oaat of XI la diffarant* 
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and the step becomes progressively steeper for Ni and Cu, 

This effect can be understood by noting that for Pe the 

state at Ep has dominantly d— character while for Ni and 

Cu the character progressively changes to p-like. This 
behaviour is a consequence of the hybridization of the 
two ^ bands. For V and Cr/ the high moment'um (p > 2 a.u) 

component (HMC) are absent while for Pe, Cu and Ni they are 

appreciable even at p = 3 a.u. and display smaller steps 
arising out of the ISmklapp-processes (Chapter l). Contribu- 
tion to HMC mainly comes from the d-like states. The non- 
occupation of the d-like states around in V and Cr 
explains 'the absence of the HMC in these metals. In Fe, Ni 
and Cu, the number of d electrons being more, the higher 
band starts getting occupied and it contributes to the HMC 
components (Pigs. 5,1 and 5,2). 

The half width of the EMD contributed by the lov/er 
band decreases from 0,85 a, u, (v) to 0. 6 a.u. (Cu) and 
the peak value p(0) increases from 0,68 (V), to 0,89 (Cu) , 
This behavioxir of P (O) signifies that the s-like states at 
become more and more spreadout in r-space as one goes 
along the series. The effect may also be due to the- decrease 
in the d bandwidth and sp-d hybridization. Although the 
hybridization effect is seen for Fe, Ni and Cu around 
p = 0.4 - 0.8 a.u. in the Et4D ciirves of -the individual bands, 
the total EMD does not show such effect significantly. 
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Flowever, in • V and Cr where only the lowest hybridized band 
contributes^ a small bulge is seen near the zone face* 

b) ^ (110) Direction ; 

This is the most interesting direction to compare the 
systematics of EMD, Firstly^ the 154D along this direction 
shows I-MC for all the metals. Secondly the higher band lies 
close''- iio the Fermi level/ thus making the U-lD sensitive to 
the structxire of Fermi surface. The HMC arise from the partial 
occupation of the higher 2^ bands. Fig, 5,7 shows how the 
position of Fermi level affects nature of EMD in V [lO], in 
Fig, 5.7/ index 1 describes the case in which Fermi level 
cuts £ ^ band twice giving rise to a miiLtiply connected 
jungle-gym surface Consisting of interconnecting arms along 
(lOO) direction and close hole pockets around K. The corres- 
ponding EliD shows remarkable structure. On the other hand if 
V 7 e take to lie above point (Case 2')/ then the IMD shows 
no such structure, 

Wakoh et al» [ 11] have shox'Tn that experimental results 
support case 1. Similar analysis can be made for Gr. The 
presence and size of the steps in Fe, Ni Cu can be explained 
in terms of the d— character at Ep as in the case of ^ direc- 
tion. 















A 


L r' 


b«adB and jP (p ) for Hi and Ou 
dloa^ fill] diraotlAa* ^Ria T»rtle#a 




case 1 
case 2 



as disouse sd in the text. 
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Along the I direction the I-JVIC component is shown up 
even for V and Cr, The contributions to these HMC come 
mainly from the U— processes due to the 'd* like character 
in the bands, 

c) ^ (ill) Direction : 

The A direction can be analysed in a similar fashion, 

A peculiar feature of the bcc metals is that, all of than 
show a sharp drop in the HVID in the 2nd zone due to 
U-process; even though the lowest A^ band is oompletely 
occupied. This happens because the contribution to the 
EMD from P to H' in the second zone comes from the U-process 
with the k lying on the PH line (or P bands), where P^ is 
being partially filled. In Ni and Cu the EMD shows a smooth 
behaviour since both the bands are completely occipied, vrhich 
is associated with the 'neck* of the PS, 

It may be added that we have also calculated the EMD 
in Rh and it shows a trend similar to the 4d (fee) metals 
Pd and Ag. A systematic study of the 4d metals cannot be 
made until the EMD in lib; Mo is calculated. Work in this 
this direction is now in progress. 

To summarise, the analysis of this chapter indicates 
that the anisotropic contribution to the EMD, P(p) , comes from 

i) the Fermi surface, and 

ii) the d-electrons. 
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The effects of the PS are dominant in the first zone. 
For the high moment-um region the anisotropies come mainly 
from the d electron wavefunctions/ and their t2g and e^. 


character 
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ChapteJL 6 


PTRECTIOWAi COmOH PROFILES AMD AmOTROPIES 
OF Pd AND Ag 


6,1 Int.roduct.ion: 


In the last two chapters a detailed st^Jdy of the 
EMD in transition metal S/ partictilarly in Pd and Ag was 


presented in terms of their respective band structure. As 
shown by Eq, (1,7), the Compton profile (or the ACPAR) does 
not measure the EMD, p (p), directly. It was noted in 
Chapter 1 that a comparison between theoretical and eKperi— 
mental results can be carried out in two ways: (i) recons- 
truct three-dimensional distribution P (p) from measured 
profiles and compare it with theoretical P (p) or (ii) construct; 
directional profiles* from theoretical p(p) and I 

compare these (or their differences) with experimental resifLts* ■ 


The former method of reconstruction of EMD has been 
suggested by Mijnarends [l ] who expanded the EMD in a series 
of lattice (Kubic) harmonics. Einploying this method, direc- 
tional profiles from single crystal specimens of Cu [l], 

Fe [2], Li [3], diamond [4], Si [ 5], Fe^O^ [ 6] , [7], Mg [8], 

Cd [8] and Ni [9] have been analysed. Recently Mueller [10] 
has developed a Fourier method to reconstruct EMD and has 
applied it to Si, 

* In the following we shall use the coninon teirm 'profile for 
the directional Compton profile '^hkl^^z^ the ACPi® curve 
^hkl^Pz^ and denote these by a common symbol 
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The above method of reconstruction of EMD is beset 
with some problems. Majumdar [ ll] has analyzed the method 
of Mijnerends and has pointed out that truncation errors can 
be serious and that unique solutions of p(p) cannot be 
obtained unless the truncation coincides with the correct 
number of terms in the expansion. For all the cases [1-9] 
analyzed, the reconstructed momentum density shows an 
oscillatory character at low momenta and only in some cases 
the tl'ieory and experiment are compared satisfactorily. 


We feel that the other approach of constructing the 
directional profiles theoretical P (p) is equally 

important as demonstrated for V [ 12] , Cr [l3,14], Fe [14], 

Ni [14] and tJb [15]. These authors have obtained by 

integration of P(p)/ as shown by Eq. (1,7) in the k-space. 
However, an integration in ]c-space dcanands a knowledge of the 
P (;^ at many p-points, thus involving considerable computa- 
tional effort. Instead we have explored another approach, 
based on an expansion of P(p) in terms of lattice harmonics. 
This method is discussed in the next section. 


A comparison between CP theory and experiment can be 
done in a more sensitive manner if v^e compare the difference 
(or anisotropy) curves " '^h'k’l* rather than the 

directional profiles This can be seen by noting that 

most of the core contributions (about 50% in transition 
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P (;^ at many p-points, thus involving considerable computa- 
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directional profiles This can be seen by noting that 

most of the core contributions (about 50% in transition 
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metals) as well as some of the experimental uncertainties 
are subtracted out in the difference curves, making them 
more sensitive to the band electron wavefunctions* Keeping 
this in mind we have presented only the difference or 
anisotropy cxirves for Pd and Ag in this Chapter. 

6*2 Spherical Harmonic Expansion of BMP ; 

i’ollowing Mijnarends [ 1] we expand p in a series 

of lattice harmonics (n^) of the appropriate symmetry 

P 

a V p 

where the P (n^) are an orthonormal set of invariant 

p 

linear combinations of spherical harmonics of order A and 

Q ^ specifies the orientations of p with respect to a suitably 
P 

chosen orthogonal coordinate system. The index v distin- 
guishes the various harmonic of the same order and symmetry 
type. For cubic lattices which we are considering, only 

harmonics v^ith a = 0, 4, 6, 8, 10, 12 (2X), 14 will 

contribute [l6] , and m = 0 (mod 4) and positive. 

It is now mor^ convenient to define two orthogonal 
coordinate systems and R^, The first system R^ = ( 5/ n ,5) 

is fixed to the crystal and provides the reference system for 
defining the lattice harmonics. The second coordinate 
system R^ r. (x, Yt is connected to the apparatus so that 
its 3-axis is along the scattering vector. The two systems 
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and can be brought into coincidence by a' rotation of 
through the Euler angles [17] ot^ g ^ Y* The angles 3 and y 
measure the polar and azimuthal angles of the z-*axis (i,e, 
the scattering vector) respectively. The angle a measures 
the rotation of the specimen around the scattering vector. 


The lattice harmonics P ( St ) are a linear combi- 

p 

nation of spherical harmonics 


P, (n ) 

Jlv 

P 


E a (n ) 

^ £ mv A ^ 


( 6 . 2 ) 


m 


and as shown by Mijnarends [l ]/ the directional profile is 
given by 

Ab V 

where 

Siv'Pz’ ' !“ <Pz/p5 P 

iPa I 

Equations (6#l) to (6* 4) constitute our scheme for the 
calculation of the directional profiles Jg a ^ 2 ^ from p(p). 
We proceed through the following steps [ 18] to achieve our 

objectives. 

i) • Employing band structure method we calculate the 
energy bands and the wavef unctions# 

ii) We then calculate the EMD, P(p.) , for several directions 

->> 

of p. 
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iii) Using the expansion (6.1) we express p(p) in terms 

of and in a set of linear equations. 

P 

iv) We then solve these sets of equations either by direct 
inversion or by using a least-square fit/ taking cart^ that 
the number of terms (1 - values) retained should be sufficient 
to ensure convergence* 

v) Knowing we obtain g, (p ) by (6.4) and then 

proceed to obtain (p^) by (6,3). 

6,3 Application to Ni and Gu ; 

To test the above procedure we first applied it to 

the EMD of Ni and Cu for which ejqjerimental data about the 

anisotropies observed from CP and ACmR measuranents were 

available. Theoretical results of EMD^ P , used for this 

worJc were taken from Singru and Mijnarends [l9]/ Singru [20] 

and from our work. Theoretical EMD, P (p), is usually 

calculated for a limited number (n, ) of directions ft ^ , 

P 

Therefore, in practice, the expansion series (6.1) is 
truncated after a finite number (n^, where n 2 < n^) of 
terms* Since the theoretical P (p) for Ni and Cu were 
available only for the three symmetry directions [lOO ], 

[110] and [ 111] / we! '-have used n^^ = 3 and n 2 = 3 for our 

analysis of Ni and Cu. In the next section, we shall 
discuss the effect of different values of n^ and n 2 on our 
results. 
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Figs, 6.1 (a) and 6.2(aO show the anisotropies for 
the CP and ACPAR curves for Ni as calculated by us with the 
method outlined in the previous section. Experimental CP 
curves measured by Eisenberger and Reed [21 ] are shown in 
Fig, 6,l(b) while Fig, 6*2 (b) shows the ACPAR curves 
reported by Kbntrym— Sznajd et al*. [9 ] , Theoretical aniso- 
tropies for the CP and ACPAR curves are shown in Figs. 6, 3 (a) 
and 6,4 while the experimental CP curves by Eisenberger 
and Reed [ 2l] and experimental ACPAR curves by Cushner et al, 

[ 23] are shown in Figs. 6,3 (b) and 6,5 respectively, 

A visual comparison of the theoretical and experimental 
anisotropies* in Figs, 6,1 to 6.5 reveals that although the 
main qualitative features of the obseirved anisotropies are 
reproduced fairly well by our theoretical curves, there is 
poor agreement in some areas. These observations are 
summarised below. 

a) In general all the curves show a poor agreement in 
the very low momentum region p = O — 0, 2 a,u. Most satisfac- 
tory agreement is seen for the difference curves "^^lO "* *^100* 
The agreement worsens as one goes from *1210 ” '^111 ^ 

"^10 0 ■■ ■^ 111 * 

b) Theoretical curves reproduce the observed 

maxima and minima with a reduced intensity of the lobes. 

* The different momentum scales used should be noted while 
comparing theori' and ^cperlment* 













Difference in arbitrary units 
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c) The agreement for Gu (CP and ACPAR) is somewhat better 
than that for Ni. 

d) Compared to the theoretical curves, experimental 
ACmR curves show a struct\.ire at low momenta. 

To understand some of these f eatiores we should note 
from Eqs. (6.3) and (6.4) that the behaviour of parti- 

cularly at low momenta, can be affected by the following two 
types of truncation errors. 

i) Truncation errors in the harmonic expansion P. (p): 

A»V 

Arbitrary truncation in Eq. (6.1) can lead to 

incorrect representation of p (p) through p^ (p) , Again 

a choice of n^^ == 3, n 2 = 3 may allow us to fit the function 

p^ (p) at specified momentum (p) values but it may not lead 

to correct Pj, (p) values at intermediate p-values. The 

step-function character in p (p) due to Fermi surface can also 

lead to incorrect P (p). 

i 

ii) Truncation errors in the integral for g (p^): 

i,v 

The calculations for p (p) were made only rpto 

P = 3,5 a.u, because of the limitations of the Hubbard- 
max 

Mijnarends method and computational problems at high 
momentum values. This would introduce a truncation error 
In (6.4) because the upper limit of integration is brought 

down from infinity to P^^ax* 
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The latter type of truncation error can be particu- 
larly more serious for CP which have significant high 
momentum components. On the other hand for ACPAR the 
former type of truncation error can cause more problems. 

We feel that the disagreement observed in Pigs, 6,1 - 6.5 
is due to the above truncation errors, although all the 
features may not be explained so easily. However the 
results for Ni and Cu have provided us vjith a valuable idea 
about the limitations of the above method, 

Aioove analysis for Ni and Cu show that it is worth- 
while to extend our above procedure to the theoretical p(p) 
of Pd and Ag, if higher harmonics are included in Eq, (6,l1>, 

In the following section we present our results and discuss 
than, 

6.4 Results for Pd and Ag and Discussion t 

Out results for the CP anisotropy for Ag losing the 

p(p) results of Chapter 4 are shown in Figs, 6,6 and 6.7, 

We have tried to make systematic study of the truncation 

errors involving P (p) by making calcifLations for different 

a 

combinations of n^ (number of directions of p for which 
^ ip) is calculated) and n 2 (number of terms retained in 
harmonic expansion (6,l)). The values of n^^ and n 2 used in 
the analysis are marked on each figure. Similar plots of the 
CP anisotropy for Pd are shown in Figs. 6.8 and 6.9. It can be 
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observed that the anisotropies for Ag are very similar to those 
of Cu (cf. Fig, 6,3) as expected from the similarities in 
their EMD (Chapter 5). On the other hand^ the anisotropies of 
Pd are not altogether similar to those of Ni, an observation 
already made about their EMD’s in Chapter 5-. Examination of 
Eigs. 6,6 and 6,7 shows that for Ag anisotropy curves change as 
^2 changed from 3 to^ 6 , The changes are most significant 
while going from n^ = 3, n 2 = 3 to n^ =12, n 2 =5. Changes in 
going from n 2 = 5 to n 2 = 6 (n^ being 12 ) are minor indicating 
that inclusion of 6 th term in the expansion has little effect 
on the anisotropy, other observations can be stmmarised as 
follows, 

a) The difference curves " *^ 100 ^ change least as 

(n^, n 2 ) changes from (3, 3) to (12,6), It may be recalled 
that this difference curve shov;ed the best agreonent (among 
the curves discussed in Sec. 6,3) between theoiy^ and esxperi- 
ment for Ni and Cu, 

b) The difference curve (J^io *” '^111^ seems to develop 
a well-defined minimum around p = 2’ a.u, as (nj^, XI 2 ) changes 
from ( 3 , 3 ) to ( 12 , 5 ) or (12,6), 

' c) The difference curve (J^qq - shows maximrsn 

changes as ^^ 2 ) is changed from (3,3) to ( 12 , 6 ), There 

seems to be a decrease in the intensity of the peaks and a 
development of oscillations, it may be pointed out that 
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fox" Ni and. Cu/ "the diffenance curves ^ had. 

sho'wn leasL satisfactory agreement with the experiment. 

d) The effect of including higher terms (n 2 ) in the 

expansion seems to broaden the pealcs in the curves. Our 
calculations (not shown in Figs, 6,6 — 6,7) further 
indicate that the differences at low momentim region 
(around p = 0) fluctuate as n 2 goes from ’7 to 11. In the 
absence of any experimental data to compare with our results# 
we feel that for n^ = 12, a reasonable choice for n 2 is 
n 2 = 5 or 6, 

A further discussion of our results for Ag would be 
made in terms of Pj, and g^ functions later in ^±s section. 

In the case of Pd we have calculated P (p) along 
7 (n^ =7) directions only. Our observations for Pd can 
be summarised as follows. A significant change in all the 
three difference curves is observed as (n^, n 2 ) is changed 
from (7,4) to (7,5). Such a significant change was not 
observed by us for Ag, It is known that for a reliable 
1 east~square*-f it as in the inversion of Eq. (6,l), we should 
have n^/2. This might explain the above behaviour.^ In 

any case the EMD in Pd is more sensitive to the truncation 
errors in the harmonic expansion . (6,1 ) . Other changes in 
Pd caused by (n^, n 2 ) going from (3, 3) to (7,4) are quali- 

tatively similar to those observed in Ag. 
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A deeper insight into the above results can be obtained 
by studying the behaviour of the functions aj(p) and g, (p ) 
in (6,l) and. (6.3) nespectivaly. 

For this discussion it is helpfioi to knox-j- the angular 
dependence of the lattice (kubic) harrnonics in the (lio) 
and (lOO) planes as shown in Fig, 6,10, The • and g^^ 
functions obtained by us for Pd and Ag are shown in Figs, 6. 11 
to 6,13^ respectively. These functions reveal the following 
important features, 

a) The functions for Pd show more structure than 

those for Ag, Functions represent the spherically 
averaged part of the EMD, P(p)/ while higher P are a 

JO 

measxjre of the anisotropic part of the EMD, In Chapter 4, 
it was already noted that the EMD in Pd is more anisotropic 
than that in Ag, This fact seems to be reflected in the 
behaviour of the respective P^^ functions. The higher (Jl ^4) 
terms in Ag show contributions around the Fermi momentvim 
Pp and around pp + G, but in Pd they have appreciable con- 
tributions at intermediate momentxam values also. This 
indicates that the anisotropies in the MD of Pd are built 
up not only from the Fermi 'sxirface anisotropy but also from 
a strong ' :k— dependence of the d— electron momentxm distribution. 
The curves also indicate that the EMD in Pd is more compli- 
cated than the EMD in Ag. Other features displayed in 
Figs. 6,11 to 6,13 can be xanderstood from the resxxLts of the 
EMD of Pd and Ag taken together with Fig. 6,10. 




{a) Pd 


functions 


(b) Ag 




go(p) in arbitrary units 



Fl«, «.13 

ft) S^pliftrieft3l3r ftw %) ay^siftiiSl.y' ariarM^ %KX> 

p^sift, i^(p), Pd m& r^t)» tw M tm Ag, 
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B} For Pd, Pg shows an appreciable contribution xinder- 

lining the need to consider higher harmonics ( ^ i 8). 

c) The function g^ represents the spherically average! 

part of the CP v/hile higher ( £ 2 ^ 4) g ^ functions give 
the anisotropic part of the CP. The plots of 9, (p ) 

Fig. 6.12 display a rather slow convergence. 

To sxmm arise, the analysis made in this Chapter has 
shown the following. 

i) Starting from a theoretical EMD and employing the 
lattice harmonic expansion metlxjd outlined above, a fair 
idea about the anisotropies of the directional profiles 

can be obtained. This information can provide a preliminary 
guideline for experimentalists in selecting the particular 
single crystal directions (h, k, l) which could be most 
interesting. Of course the best vray to compare theory with 
experiment is to carry out a two-dimensional integration of 
p (p) in momentum space, tfowever the lattice harmonic 
expansion method can offer a faster but approximate way to j 
obtain a rough idea of the anisotropies expected or to | 

compare theory with experiment, [ 

: i 

ii) » We feel that the main utility of the present method 
is that an examination of the and g£ plots can be ! 
instructive for carrying out an 'unfolding* [ l ] of measured 
CP (or ACffi-R) data for reconstruction of P (p) * Thus our [ 
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graphs (Figs, 6.11 to 6. 13") suggest that for Ag and Pd one 

shoiiLd measiare CP along enough number of directions so that 

at least 6 terms (upto 1 ^ = 12) can be used in the hartKjnic 

max 

expansion. 

A better xanderstanding of the usefulness and applica- 
tions of the present method can be obtained once experimental 
CP data for Pd and Ag becomes available. We feel ’that a 
folding and unfolding analysis by the lattice harmonics 
expansion method on same single crystal metals should be 
carried out to gain -a better insight into the method. It ma^ 
be worthwhile to attempt a similar analysis employing the > 
Fourier method [ 10] as well. i 

Our analysis has also brought out the limitations of | 

the lattice harmonic expansion method. Firstly# the disconti-j 

■ 

nuities in the P (p) cannot be handled satisfactorily by I 

I 

this method. Secondly the choice of (n^# to be used is [ 

not unique unless n^ is quite large as compared to n 2 (e. g.i 
n- s 40# n,, * 10)# thus necessitating large computational { 

12 j 

efforts. Thirdly the convergence of g^ function is slow i 
and a certain anount of arbitrariness always ranains about | 

■ i 

the truncation in (6#l), ; 

I 

i 

" j 

■ ■ '■ f 

■ , ■ ■ f 
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